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1 Introduction 



The d = 2 minimal conformal quantum field theory whose central charge 
is equal to c = 7/10, known as Tricritical Ising Model (TIM), describes the 
scaling region near the tricritical point (i.e. the end point of a line where 
three phases coexist simultaneously) of the universality class corresponding 
to the Landau-Ginzburg theory with an interaction 0^ In particular 
it can be thought of as a continuum realization of the Ising model with 
annealed vacancies, see and it describes adsorbed hehum on krypton- 
plated graphite, 0, (for a general reference see 0]). 

As is well known in the d = 2 minimal conformal theories there exists a 
finite number of primary operators In TIM they can be divided in two 
classes according to their parity with respect to a natural Z2 symmetry: the 
energy operators are even and the spin operators are odd. At the critical 
point there is a spontaneous breaking of the Z2 symmetry. TIM displays 
other symmetries apart from the conformal one: super conformal symmetry 
1^, 13 and those based on the coset formulation (su(2) and 67 algebras, see 

i)- 

In this article we study the two point correlation functions of the energy 
operator and of the spin operator in proximity of the tricritical point. More 
specifically we consider the model obtained by adding to the conformal action 
the most relevant energy operator of the critical theory, e (of scale dimension 
1/5, see Section ||): 

A = AcFT + A / d^xe{x). (1.1) 



We also restrict ourselves to the high temperature phase of the model, that 
corresponds to positive values of the coupling A. 

The considered model is integrable and its spectrum and S matrix are 
known, |]ll| (see also IQ)- In the high temperature phase the massive excita- 



tions are given by seven elementary particles. The mass of the fundamental 
particle is given by: 

m = Cf\^'\ (1.2) 



where Ci? is WM: 




4 ^^r(|)r(f; 



5/18 



)'r(|; 

3.745372836243954223- •• (1.3) 



The masses of the other particles can be found in [|1^ with a more detailed 



analysis of the model. In the same paper the authors were able to find 
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the form factors of the model (see [T^ for general references and |jT6| for 



recent developments) and to calculate the energy-energy correlation function 
including the excited states given in table 3 of that work. The representation 
of the correlation function in terms of form factors gives rise to a long distance 
expansion (the small parameter being e~^^) that is well appropriate in the 
region down to mr ~ 1, but is reasonably less quickly convergent in the 
mr << 1 region. 

In this paper we study the energy energy and spin spin correlators by use 
of a perturbative expansion around the conformal theory: in this case the 
small parameter is mr and consequently our analysis is complementary to 
the form factor approach. 

The perturbative calculations are performed by using the Operator Prod- 
uct Expansion (OPE) approach in order to overcome the infrared divergences 
problem typical of the relevant perturbations of massless theories. The idea 
of this approach came out very early in the general context of quantum field 
theories |17| (see for similar ideas in QCD), was first introduced in the 



context of perturbed conformal theories by Al. B. Zamolodchikov in (see 



also pO| for an analog independent proposal and for an early application 
of Ref. [0), and was developed at all orders in a very general context in p2| . 

In the OPE approach the terms non analytic in A, whose very existence is 
related to infrared divergences of naive perturbation theory, are segregated 
in the Vacuum Expectation Values of the operators present in the OPE 
expansion of correlators 

< $„(r)<|.,(0) >A~ E Cl,{r, A) < $^(0) >a (1.4) 

c 

($a are deformations of the conformal theory operators, the suffix A refers 
to the complete theory correlators; the dependence on A of the Wilson coef- 
ficients will be omitted in the following). The resulting perturbation theory 
for the Wilson coefficient C^f^{r, A) turns out to be IR finite at all orders. The 
required operators VEV are a non perturbative additional information that 
should be considered as an external input in this context. 

In the specific example that we consider, the mean value of the energy 
operator e will be determined by the Thermodynamic Bethe Ansatz, |T2 



and the mean value of the sub-leading energy operator, required for our 
purposes, will be estimated to a reasonable accuracy by using two exact sum 
rules that express the scale dimension of the energy and the central charge of 
the unperturbed conformal theory in terms of moments of the energy energy 



correlators |T0| (in the long distance region we use the results of |]T3[ for 
the correlator). 
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Our final result is an estimate of the energy energy and spin spin corre- 



lators up to 



mr] 



The prediction for the energy energy operator will then 



be compared with those of [jT3| in Figure 1. 

The plan of the paper is as follows: in Section |^ we collect the required 
knowledge on the TIM and introduce the first order perturbative formulas 
for the Wilson coefficients; in Section |^ the final 0((mr)^) expression for the 
correlators are reported; in Section ^ the missing VEV of sub-leading energy 
is estimated by use of exact sum rules and long distance results. Details of 
calculations are reported in [Appendix B| while [Appendix A| reports a very 
general theorem that turned out to be useful in calculations. 



2 Tricritical Ising model 

The Tricritical Ising Model (TIM) is the minimal CFT with central charge 

). At the critical point we have the 



10 V"- ^ t' notation of 
following operators: 



•the identity $11 = /, A/ = 0, 
»the magnetization field $22 = o", ^o- = Jj, 
»the sub-leading magnetization operator $12 
»the energy density $21 = e, = j^, 
» the vacancy operator or sub-leading energy operator $31 = /, A; = |, 



^1 j^g , 



the irrelevant field $ 
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s", A, 



Using the fusion rules and knowing the Wilson coefficients we can con- 
struct the algebra for TIM: 
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The corresponding Wilson coefficients can be found in Table 4 of p4 
(their normalization: < >= 5ij will be used everywhere here.) 

Let us write here the adimensional Wilson coefficient C'^ (C = Clrl*^™*"): 



CI 



2 r(4/5)P(2/5) 

3\ r(i/5)r3(3/5) 
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= 0.6103020742648446931 ■ 



(2.1) 



and 

= (2.2) 

We consider in this paper the following perturbation: 

A = AcFT + A y" d^xe{x) (2.3) 

where A ~ T — Tc > and e is the most relevant energy operator of TIM, 
with total dimension 1/5. 

The idea of the OPE approach (see |T^, [1^, [2^]) is to confine contributions 



non analytic in A inside the VEV's and to build up an IR finite perturba- 
tion theory (Taylor expansion around A = 0) for the Wilson coefficients. 
Assuming the regularity of the Wilson coefficients in terms of the renormal- 
ized coupling A, the asymptotic convergence of OPE and the validity of the 
Renormalized Action Principle, that in our case looks like 

- ^ <[■■■] >^= /< N •■■■] ■e:=e-<e>x (2.4) 



([■••] means renormalization, when needed), an IR finite representation for 
the generic n*^ derivative of C^^ with respect to the coupling, evaluated at 
A = 0, has been given involving integrals of conformal correlators. See 



23 1 for an application to the critical Ising model perturbed by a magnetic 



field. 

In our case we will approximate the correlators by use of the following 
first order expressions (all derivatives below are intended to be evaluated at 
A = 0): 

<ee>x = Ci + C',,<l>x+\dxCl<e>x+Oi\^) (2.5) 
<aa>x = Cl„ + XdxCl^ + Cl„<s>x + 

XdxCl^ < e >x +Cl„ < I >x +0(A^/3). (2.6) 

At first order we derive the following nontrivial relations, (see also |[T9| , [20[ 
for first order formulas): 

- dxCl = fdh < s{z){air)a{0) - CtAO)) > (2-7) 

-dxC^^^ < 5(oo)£(0) >= 

d'z < e{oo)e{z){a{r)a{0) - - > (2.8) 
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-dxCI, < e{oo)e{0) >= 

d'z < e{^)e{z){e{r)e{0) - C], - 0^(0)) > (2.9) 



where the prime means any (rotation invariant) IR regularization of the in- 
tegrals. We observe that the derivative of Wilson coefficients is obtained by 
adding up a " naive perturbative term" (insertion of the interaction J e in the 
correlator) and some "IR counterterms" generated by the theory itself in a 
natural way: both contributions are IR divergent but their sum is finite and 
gives the wanted result. This structure persists in the general all order for- 
mulas. The explicit form of the correlators and the details of the calculation 
of the IR finite part are reported in [Appendix B. 



We can parameterize the VEV of e and / as 

<e>x = A,X'/^ (2.10) 
</>A = AiX^/\ (2.11) 

Notice that potential adimensional logarithms are absent in the VEV's be- 
cause no renormalization effects arises in this case (an operator of dimension 
X can mix in this theory only with operators of dimension x—^k for positive 
k). 

The constant can be fixed by use of results of [|I^ (Thermodynamic 
Bethe Ansatz) 

A - 5 sil) ( An^ninif Y' ( 2r(i) y 

' " 36 5(1)5(1) [ r(i)^r(|) ) [nim) ^'-''^ 

= -1.468395424027621489- •• (2.13) 

where S{x) = sinvrx (see e.g. for more details). The quest for an estimate 
of Ai is the goal of Section H 



3 Results 

We report in this Section the results for the approximate energy energy 
and spin spin correlators, Eqs.(|2.5|-P^). More details on the calculations can 



be found in Appendix B 



It is convenient to define adimensional quantities C = C|r|*^™*-^, Fa 

< AB >h\r\^^^^^ and introduce the scaling variable r = mr, where m is 
the fundamental mass of the model, defined in Eq.( |1.2D . 
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We have: 



1/5 _ / X 6/5 



= l + A,Ci(f] +0(ri«/^) (3.1) 



1 + A,CIA—^ +AiCL[-^^ 



+9aC1,(— j j +0(r3) (3.2) 



By use of results in [Appendix B we obtain: 



s{\)'s{^'s{^s{^) 26/5r(i)^r(|) 

= -1.008826378183031- ■■ (3.3) 

^(|) r(|)^ lvr(|)r(^^' 

= -0.223579589302507- ■■ (3.4) 

^ .^ j(ifmfs{m^fs({f (-g(t) + sci) + s{f) + sci) - sci)) 

s{lfs{lfs{l)sCi)sCi) 2W5r(i)^ 

= -0.266530272188332- ■• (3.5) 

All other quantities have been defined previously; for a numerical value of Ai 
see Eq.(^). 



4 Sum rules 

In this section we will consider two known sum rules satisfied by the exact 
complete theory. 

The sum rule more sensitive to the short distance behavior of the corre- 
lators is the one concerning the scale dimension of operators. If we consider, 
as is the perturbation A / e of a conformal field theory (assumed to 

draw the system towards a massive theory) and an operator X with scale 
dimension Ax, assumed to have a VEV < X >x= AxX^^^^^~^^\ then, by 
taking the derivative with respect to A of the VEV and by using the Action 
Principle hypothesis, Eq. (|2.4|) , it is easy to obtain 

= -A / d^x < [X{x)e{0)] >xe/ <X>x. (4.1) 
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This exact sum rule has been derived first, in a more generahzed form, in 

m. 



By choosing X = e and using the known relation between the trace of 
energy momentum tensor and the perturbing field e (holding at least when 

A, < 1/2), 

e = 47rA(l - A,)e (4.2) 
we have the equivalent relation: 

d'x<QQ>,^=--^^m'. (4.3) 



The second sum rule we consider is given by 



c= ^ = ^ I d''x\x\'' <ee>xc (4.4) 
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and has been obtained in the context of the C-theorem 

The short distance behavior of the connected 00 correlator can be ex- 
pressed in terms of the scaling function of Eg. (pTTP as 



(Notice that we subtracted < e >\ to get the connected correlator). 
The long distance behavior is given by: 

< 00 >xjm^ ~ -{f^Ko{c2r) + f^Koic^r)) (4.6) 

TT 

where cs = 2 cos(57r/18), C4 = 2cos(7r/18) and /a = 0.9604936853, U = 
-0.4500141924, HTl. 



We now proceed in the same way as done in [^], where good results were 
obtained: we split the integrals in the sum rules at some point mr = A, 
evaluating mr < A region with the short distance expression Eqs. p.l| ),( [4.5D 
and mr > A with Eq. ([4.6|) . In particular the integral of the long distance 
correlator is given by 

/ < 00 >Ac rfV/m^ = 2A(/|/c2iri(c2A) + fJc,K,ic,A)) (4.7) 

J\mr\>A 

while for the c sum rule the integral of the long distance contribution is given 
by 

/ < 00 >A c r'd'x = 2(/|/c^i?(c2A) + fJctKicA)) (4.8) 

J\'mr\>A 
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where 

K{x) = {Ax + x^)Ki{x) + 2x^Kq{x). (4.9) 

By imposing the sum rule Eg. ( [4 .31) we can recover the unknown constant 
Ai as a function of A. Following a good approximation for Ai can be 
obtained by looking for a minimal sensitivity point A* of Ai. We obtain in 
this way 

Ai = 3.78 (4.10) 

where an error of 5 — 10% is extimated by using different variants of the 
method (imposing validity of Eq. ([4.4|) or approximate validity of both). 



5 Conclusions 

By use of results of Sections we can give a prediction on the short 
distance behavior of the energy energy (or BG) and spin spin correlators of 
the TIM model near criticality (most relevant energy perturbation). As a 
by product of our work we get an estimate for the subleading energy VEV, 

Eqs.(^)-i;Tg) 

The most striking feature of OPE approach used in this paper, is the fact 
that non analyticities in A, present in the short distance behavior of corre- 
lators, are naturally recovered from the VEV and manifest in non negligible 
fractionary power corrections (in A) in Eqs.( p.lD - ( p. 21) . 

In Figure 1 we compare our predictions for the GG correlator as a function 
of mr with the corresponding long distance prediction of [0 . The agreement 



between the two approaches is good in the intermediate 0.5 < mr < 1.5 
region, where corresponding "small" parameters are actually of order 1! See 
M for similar situations. Moreover the lower mr region shows a 



reasonable evidence of convergence of the long distance expansions towards 
our results. We remark in particular, referring also to an similar comparison 



made by the authors for the critical Ising Model in magnetic field [23, fig.2], 
that the weaker the singularity in r = of a correlator is, the fastest is the 
convergence of the long distance approximants. We think that, in any case, 
the two approaches are always complementary. 

Finally we want to stress that the OPE approach gives an all order IR 
finite perturbation theory that is very general (mild and reasonable hypothe- 
ses are behind it) and do not uses integrability essentially: the computation 
of Wilson coefficients can thus reasonably applied to many others physically 
interesting models in statistical mechanics. The main problem of reaching 
higher orders is related to the complexity of the required integrals. We think 
nevertheless that second order expressions when necessary, could be reached. 
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at least numerically. The ignorance of the expressions of the VEV will result 
at each finite order in a finite set of universal constant parameterising all the 
correlators. How to reach a non-perturbative information on those constants, 
apart from what already known from Thermodynamic Bethe Ansatz, is still 
an open problem. 
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during a visit of the authors to CERN, whose warm hospitality is acknowl- 
edged. The work of R.G. is supported by a TMR EC grant, contract N° 
ERB-FMBI-CT-95.0130. R.G. also thanks the INFN group of Genova for 
the kind hospitality. 



Appendix A A useful theorem 



We give here a (slight generalisation of) a useful theorem proved in p8| , 
(see also for related work). 



Let us consider a double integral on the plane of the form 



/ 



N 

d'w J2 fMQM^n (A.l) 

a,/3=l 



where {/q(w)}^=i {//3('W^*)}^=i are two sets of linear independent functions, 
Qaf3 is a complex, constant, matrix and * denotes the complex conjugate. 

We will assume that fp{w*) and {f[^{w))* have the same monodromy 
structure. In particular we assume firstly that the two sets of functions fct{w) 
and gpiw) = {ff3{w*))* have the same (algebraic) branch points, {wk}^=Q, 
such 

= \Wo\ < \Wi\ < ■ ■ ■ < \Wm\ < \Wm+l\ = OO (A. 2) 

and that are analytic elsewhere. (The relative cuts in the w plane are taken 
to stay along the same line 7 connecting the w^s and chosen such that \w\ 
is increasing along it. Cuts in w* plane are obtained by * mapping the 
cut w plane.) In the following we will require also that for each choice of 
k = 0, ■ ■ ■ ,m + 1, in the limit w Wk, 

\{W - Wk)fa{w)\ 

\{w-Wk)Mw*)\ ^ 0, (A.3) 
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uniformly on the cut plane. This condition guarantees the convergence of 
the integral Eg. (|A. 1| ) and the validity of some steps of the proof. 

Secondly we suppose that / and g have the same related monodromy 
matrices M^, defined by the requirement 



f{w.) = M,f{w+). 



(A.4) 



Mfc above is the same constant matrix for all points w staying in the region 
of 7 between Wk-i and Wk and + (_) means limit from above (below) the cut 
(similar definition for g{w)). In particular w± are connected by a path in the 
cut plane enclosing all Wj with j < k. 

We assume finally that the matrix Q is invariant by the monodromy group 
action: 

yk Q = MiQMl. (A.5) 

All requirements together enforce the single valuedness of the integrand 
and the stronger relation 



f{w'_)Qf{w*^) = f{w'^)Qf{w^*) 



(A.6) 



for any on the region of 7 between \wk\ and (the same for 

both). 



It follows then, after a simple generalisation of the proof of based 
on Stokes' theorem (and using strongly Eq.( [A.6| )), that we can express / in 
terms of one dimensional integrals as 



^ fc=l 



where, * means transposition. 



k+l) 



1 - M,)-' 



Q 



Cfe 



dwfiw) 
dw*f{w*) 



X, 



(fc) 



a/3 



(A.7) 



(A.8) 
(A.9) 



and Ck {Ck) are counter-clockwise (clockwise) circumferences starting at Wk+ 
{wk\) and ending at Wk- {wk*J), enclosing all the Wj {w*) of lower modulus. 
We refer to for more details (see also for related integrals). 
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Appendix B More details on computations 



The expression for the correlators in Eqs.( p.7[ - p.9|) can be obtained from 
the results of for the second order correlators of minimal conformal the- 



ories. In particular the general structure is 

< 0Mi(^l)0M2(^2)0M3(^3)0/.4(^4) >0C 

n |%|'^- {S{a + b + c)S{b)\h{a, b,c-r])\' + S{a)S{c)\l2{a, b, c;r])\' 

i<j 

where the ordered multi-index /i = [n, m) labels operators and fixes their 
dimensionality (in our case e corresponds to /i = (2, 1) while a to (2, 2)), and 

S{x) = sin(7rx) 

^12^34 



V 



2^13^24 



w , x T(-a-b-c-l)T(b + l) 

h{a,o,c;r]) = — 2Fi{-c, -a - b - c - 1; -a - c;r]) 

1 {—a — c) 

l2ia,b,c;v) = + ^)^(" + ^) a + 1; a + c + 2; r/). (B.l) 

1 (a + c + 2) 

The relation between the multi-indices fi and exponents and parameters 
a, b, c can be found in |2^ and will not be reported here. We denote with 



pFg{ai, ■■■,ap;bi,,---,bq;z) = E^o fc^tlf^q^t^^ (generalised) hyperge- 
ometric function (where (a)^ = r(a + k)/T{a)). 

The integrals we are interested in are therefore of the general form: 

z{a,b,c,d,e) = y^v^-^krii -■zri^-^zri^pii -^r' 

= S{a + c)-^ J d^r]\7]\^''\l-r]\^'' X 

(Sia + b + c)Sib)\h{a, b, c, r])\^ + 5(a)5(c)|/2(a, b, c, r])\''). 



The second line can be obtained by use of results of We evaluated 

the second part of the integral by use of the theorem reported in Appendix 



~M After the completion of the computation we realized that an equivalent 
integral could be found in Appendix D of and we checked that the results 



are the same (modulo some non trivial algebraic relations satisfied by the 
generalised hypergeometric functions, |^). We report our expression that, 
by construction, involves a fewer number of boundary integrals (generalised 
hypergeometrics) to be computed. 
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Applying the theorem of [Appendix A| to our needs, it is easy to reahze 
that in Eq.( |A.7|) there is only a contribution coming from J(i). If we pose 

j(i) _ 2i(e^"^5(rf) Ji, e^"('^+'=+'^)5(a + c + rf) J2) (B.2) 

and compute the inverse of the monodromy matrix terms, it follows after 
some algebra that: 

Z = zuJ'^ + Z22J2 + Z12J1J2 (B.3) 

Ji = f\'^{l- zyii{a,b,c,z) 
Jo 

J2= l\\l-zYl2{a,h,c,z) 
Jo 

where S~^{x) = 1/S{x) and 

= B{b+1, -a-b-c-l)B{d+l, e+l)3F2(-c, -a-b-c-1, d+1; -a-c, 2+rf+e; 1) 
J2 = B{a+1, c+l)B{2+a+c+d, e+l)3F2(-6, a+1, 2+a+c+d; a+c+2, 3+a+c+d+e 



zu = -^S{a + cy^S'^{c + d + e)S'^{a + b + c + d + e)x 
S{b)S{a + b + c)S{d) X 

{S{b - c - d) - S{b + c - d) + S{b + c + d) - S{2a + b + c + d) 
-S{b + c + d + 2e) + S{2a + b + 3c + d + 2e)) (B.4) 

Z22 = ^S{a + cy^S'^{c + d + e)S~'^{a + b + c + d + e)x 
S{a)S{c)S{a + c + d) x 

{S{a + b - d) + S{a - b + d) - S{a + b + d) + S{a + b + 2c + d) 
-S{a -b-d-2e)-S{a + b + 2c + d + 2e)) (B.5) 

Z12 = 2S{a + cy^S'^{c + d + e)S'^{a + b + c + d + e)x 

S{a)S{b)S{c)S{a + b + c)S{d)S{a + c + d) (B.6) 

(We assumed the parameters to be real as it is the case of interest). 

With the knowledge of Z{a, b, c, d, e) the IR finite part of the regularized 
"naive" perturbative integrals can be extracted by use of Mellin transform 
techniques without computing all IR counterterm integrals, as explained in 
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detail in Section 2-3. Using some known relations of 3F2 hypergeometric 
functions of argument z = 1, ||3^, and fixing the normalisation of correlators 
by imposing clustering properties (see [^), Eq. ( WM can be recovered. 



Figure Captions: 

Fig 1: Different estimates of < 00 >a c correlator plotted vs mr: contin- 
uous curve is the short distance approach (this paper); dashed curve refers 
to first two states while dots correspond to the first 11 states in form factor 
approach, [[T^ . 
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